Low- and high-frequency noise from coherent two-level systems 
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Recent experiments indicate a connection between the low- and high-frequency noise affecting 
superconducting quantum systems. We explore the possibilities that both noises can be produced by 
one ensemble of microscopic modes, made up, e.g., by sufficiently coherent two- level systems (TLS). 
This implies a relation between the noise power in different frequency domains, which depends on 
the distribution of the parameters of the TLSs. We show that a distribution, natural for tunneling 
TLSs, with a log-uniform distribution in the tunnel splitting and linear distribution in the bias, 
accounts for experimental observations. 
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Recent activities and progress with quantum informa- 
tion systems rely on the control of decoherence processes 
and at the same time provide novel tools to study their 
mechanisms. Experiments with superconducting qubits 
revealed the presence of spurious quantum two-level sys- 
tems pj with strong effects on the high-frequency (~10 
GHz) qubit dynamics. Other experiments [2J suggested 
a connection between the strengths of the Ohmic high- 
frequency noise, responsible for the relaxation of the 
qubit (Xi decay), and the low- frequency 1// noise, which 
dominates the dephasing (T2 decay). The noise power 
spectra, extrapolated from the low- and high-frequency 
sides, cross at u of order T. This is also compatible 
with the T 2 dependence of the low-frequency part, ob- 
served earlier for the 1// noise in Josephson devices |j,|4j . 
Much clearer evidence for the T 2 behavior was obtained 
recently |5j,|6|- 

In this letter we point out that a set of coherent two- 
level systems (or, in fact, arbitrary quantum systems 
with discrete spectrum) produces both high- and low- 
frequency noise with strengths that are naturally related. 
We show that for a realistic distribution of parameters 
tunnel TLSs (TTLS) produce noise with experimentally 
detected features: the 1/f behavior at low frequencies, 
the Ohmic (oc to) high-frequency noise, and the T 2 tem- 
perature dependence of the integrated weight of the low- 
frequency noise. This implies that the 1/f and Ohmic 
asymptotes cross at w ~ T as was indeed observed in 
Ref. [2j at one value of T. The distribution is log-uniform 
in the tunnel splitting and linear in the bias. Microscopi- 
cally, this distribution may describe double traps or "An- 
dreev fluctuators" considered recently by Faoro et al. Q 
in their study of the relaxation (7\ decay) of Josephson 
qubits due to the high-frequency noise. Our results are 
obtained for environments with a large number of TLSs 
which are weakly coupled to the qubit. A strong coupling 
between a TLS and a qubit can lead to resonances 0, • 

Ensembles of TLSs were discussed extensively in the 
literature. On one hand, they produce a natural model of 
1/f noise, as a result of incoherent random transitions || , 
and there is substantial experimental evidence that the 



low- frequency 1/f noise in single-electron devices may 
be produced by TLSs [9|, |lfj. In solid-state qubits, 
e.g., Josephson qubits, the pure dephasing is dominated 
by this noise [ill H^ . Ilflj . On the other hand, ensem- 
bles of coherent TTLSs were sug gest ed to explain low- 
temperature properties of glasses [lj, [l5| ■ Both "trans- 
verse" and "longitudinal" couplings, defined below, were 
discussed in relation to various physical phenomena. A 
transverse coupling of phonons or electrons to the TLSs 
is responsible for the absorption and emission of energy. 
It was invoked in the discussions of, e.g., the phonon at- 
tenuation |l6| and of the low-temperature dephasing in 
disordered metals On the other hand a longitudinal 
coupling was found to be resp onsible, e.g., for the con- 
ductance fluctuations 0, 0, H3] ■ We suggest that in 
nanocircuits, e.g., solid-state qubits, both types of cou- 
plings play an important part and produce noise with 
related properties in various frequency ranges. 

As a model we consider a set of coherent two-level 
systems described by the Pauli matrices cr p j, where 
p = x,y,z and j is the index of a particular TLS. We 
write the Hamiltonian of the set in the basis defined by 
their contributions to the relevant fluctuating quantity 
(cf. Eq. © below); 
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Here, in the language of TTLSs, Sj are the bias energies 
and Aj are the tunnel amplitudes between two states. 
Each TLS with label j is subject to dissipation due to its 
own bath with Hamiltonian fldiss,j- We do not specify 
Hdiss.j, but only assume that it produces the usual relax- 
ation (Ti) and dephasing (T2) processes. We assume that 
all the TLSs are under-damped, with Tij = T-j - - <C Ej 



and To i — T, 



2,3 = J 2j < E 3- 



hi 

Hero Ei = Je 2 + A 2 is the 



energy splitting. 

Each TLS in the ensemble contributes to fluctuations 
of a physical quantity X, e.g., the gate charge, which 
affect an experimentally accessible system and thus may 
be detected. A qubit may serve as a convenient noise 
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detector 0,H2|. E.g., in the recent experiment of Ref. 
qubits were used to investigate the properties of their 
environment. We choose 



X 



(2) 



where Vj are the coupling constants and cr 2 J = ±1 cor- 
respond to the two states differing, e.g., by the value of 
the dipole moment. The interaction of the qubit with the 
TLSs is often described via a linear in X coupling to a 
variable O qu bit of the qubit, i.e., H int (X) oc AO qu bit- 

Our goal in the following is to investigate the noise 
properties of X, that is we need to evaluate the (unsym- 
metrized) correlator 



C x {u) = dt{(X(t)X(0))-(X) 2 } 



(3) 



For independent TLSs the noise is a sum of individual 



contributions, Cx — Tlj v ' 2 j Cj, where 



/ dt{{a z>j (t)a z ,j(0))-{a z ,j) 2 } 
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To obtain Cj we first transform to the eigenbasis of the 
TLS. This gives 



#tls = E |~2 EjPzJ + H &"j} 



and 



X = 



E 

3 



Vj (cosOj p Ztj - sin0j p Xjj ) , 



(5) 



(6) 



where tan 9j = Aj / ej . The first term of |JB} produces the 
longitudinal coupling (mentioned above) while the sec- 
ond term produces the transverse one. Using the Bloch- 
Redfield theory |U 0] we find readily 
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In thermal equilibrium (p z j) = tanh(-Ej-/2T). The first 
term, due to the longitudinal part of the coupling, de- 
scribes random telegraph noise of a thermally excited 
TLS. We have assumed Tij <C T, so that this term is 
symmetric (classical). The second term is due to the 
transverse coupling and describes absorption by the TLS, 
while the third term describes the transitions of the TLS 
with emission. We observe that TLSs with Ej ^> T 
contribute to C'x only at (positive) lo = Ej. Indeed 



their contribution at lo — is suppressed by the ther- 
mal factor 1 - (p z j) 2 = 1 - tanh 2 (£' J /2T). Also the 
negative frequency (emission) contribution at lo = —Ej 
is suppressed. These high-energy TLSs remain always in 
their ground state. Only the TLSs with Ej < T are 
thermally excited, performing real random transitions 
between their two eigenstates, and contribute both at 
lo = ±Ej and at u> = 0. Such a multi-peaked structure 
of Cj{uo) was rec ently discussed in various contexts, e.g., 
in Refs. [H, 0, HJ. Note that the separation of the 
terms in Eq. into low- and high-frequency noise is 
meaningful provided the typical width of the low-w 
Lorentzians is lower than the high frequencies of interest, 
which are defined, e.g., by the qubit's level splitting or 
temperature. 

For a dense distribution of the parameters e, A, and v 
we can evaluate the low- and high-frequency noise. For 
positive high frequencies, lo ^> T, we obtain 
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N J dedAdv P(e, A, v) v 2 sin 2 9 ■ 2tt5(lo - E) 
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where N is the number of fluctuators, P(e, A, v) is the 
distribution function normalized to 1, E = Vc+A 2 , 
and tan# = A/e. Without loss of generality we take 
e > and A > 0. At negative high frequencies (lo < 
and \oj\ > T) the correlator Cx(w) is exponentially 
suppressed. 

On the other hand, the total weight of the low- 
frequency (up to u > Ti) noise follows from the first 
term of {7J. Each Lorentzian contributes 1. Thus we 
obtain 
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Here we could disregard the contribution of the last two 
terms for Ej ~ Equations JSJ and © provide the 

general framework for further discussion. 

Next we investigate possible distributions for the pa- 
rameters e, A, and v. We consider a log-uniform distri- 
bution of tunnel splittings A, with density Pa (A) oc 1/A 
in a range [A m i n ,A max ]. This distribution is natural 
for TTLSs as A is an exponential function of an al- 
most uniformly distributed parameter, e.g., tunnel bar- 
rier height |l5| . It is also well known to provide for 



3 



the 1/f behavior of the low-frequency noise [8j: the re- 
laxation rates are, then, also distributed log-uniformly, 
Pn(ri) oc 1/I\, and a sum of many Lorentzians of width 
ri centered at oj — adds up to the 1/f noise. We fur- 
ther assume that the distribution of v is uncorrelated 
with e and A. 

First, we assume that the temperature is lower than 
A max , T < A max . For the high-frequency part, T < oj < 
A max , taking the integral over A in Eq. ©, we find that 
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This is consistent with an Ohmic behavior Cx oc oj only 
for the linear density P E (e) oc e. 

Importantly, this distribution P(e, A) oc e/A leads at 
the same time to the T 2 ln(T/A m j„) behavior of the low- 
frequency weight ©, consistent with experimental obser- 
vations 0, 0, 0, 0fTlf the low-frequency noise has a 1// 
dependence, the two parts of the spectrum would cross 
around to ~ T . 

A remark is in order concerning this crossing. It is not 
guaranteed that the spectrum has a 1// dependence up 
to oj ~ T. Rather the high-frequency cutoff of the low- 
frequency 1/f noise is given by the maximum relaxation 
rate of the TLSs, ri jmax <C T, as we assumed. Then 
the extrapolations of the low-frequency 1/f and high- 
frequency Ohmic spectra cross at this oj ~ T. 

To fix the coefficients, we introduce the normalization 
constant A, so that P(e, A) = Ae/A. Then, at high 
positive frequencies, T < oj < A max , we obtain 



C x {u) ~ tt(v 2 )NAoj 



(11) 



For the total weight of the low-frequency noise we obtain 
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low freq. 



Thus we obtain a numerical factor which determines pre- 
cisely the point of crossing of the two spectra. 

In the opposite limit, T ^> A max , the high/frequency 
noise depends on the detailed shape of the cutoff of 
Pa (A) at A max . As an example, for a hard cutoff the 
Ohmic spectral density implies that P e oc e 3 , and the low- 
frequency weight scales with T 4 . For a 1/f low- frequency 
behavior, the spectra would cross at oj ~ T 2 /A max ^> T, 
in disagreement with the result of Ref. . 

Interestingly, the linear oj dependence at high frequen- 
cies and the T 2 dependence of the low-frequency noise 
can be obtained from a whole class of distributions, e.g., 
for P(e,A) = /(e/A), with arbitrary, not too divergent 
(as a function of ff) function /. Presented as a function 
of energy E and angle 8, it becomes (we have used the 
Jacobian dedA — > EdEdO) 



P{E, 6) = Ef (cote) 



(13) 



i.e., it is linear in E. This linearity ensures both the linear 
oj dependence at high frequencies and the T 2 dependence 
of the integrated weight of the low-frequency noise. 

In particular one can take P(e, A) oc (e/A) s with any 
exponent s satisfying — 1 < s < 1, and with both A max 
and e max higher than the relevant frequency oj. This in- 
cludes a uniform distribution of both A and e at all rele- 
vant energies, i.e., s = 0. For ensembles with — 1 < s < 1, 
including the uniform distribution with s = 0, high- and 
low- frequency noise is created by the same fluctuators. 
On the other hand, s = 1 is the limiting case in which 
the low-frequency noise is dominated by the fluctuators 
with 9 « 1, while the high-frequency noise by all other 
fluctuators. Yet, even in this case the strengths of the 
high- and low-frequency parts of the spectrum are re- 
lated. 

For the uniform distribution (s = 0) we obtain different 
numerical coefficients. We introduce an experimentally 
accessible constant a, such that Cx(oj ^S> T) — 2naoj. 
Then, for s = 1 (P oc e/A) we obtain from Eqs. (jl 1|> 
and H12(l that the total weight of the low-frequency noise 
is given by J u % C x (oj) « 8 ln(2)aP 2 ln(T/A min ). On 
the other hand, for the uniform distribution, s — 0, we 
obtain J l { % C x {oo) « 41n(2)aT 2 . 

We would like to emphasize that the relation between 
low- and high-frequency noise is more general, i.e., it is 
not unique to an ensemble of two-level systems. Consider 
an ensemble of many-level systems with levels \n) and 
energies E n such that the coupling is via an observable 
which has both transverse and longitudinal components. 
By a transverse component we mean the part constructed 
with operators |n)(m|, where n ^ m, while the longitu- 
dinal component is built from the projectors |n)(n|. If 
the system is under-damped, that is, if the absorption 
and emission lines are well defined, the correlator of such 
an observable will have (Lorentzian-like) contributions 
at oj — E n — E m as well as at oj — 0. As an example 
we could consider an ensemble of an-harmonic oscillators 
with X — J2j v j x ji where Xj 0,r G the oscillator's coor- 
dinates. Due to the anharmonicity Xj acquires a lon- 
gitudinal component, in addition to the usual transverse 
one. Thus a relation between the low- and high-frequency 
noise would emerge naturally with details depending on 
the ensemble statistics. 

It is useful to relate our phenomenological results to 
the recent work of Faoro et al. , where they considered 
physical models of the fluctuators, which could couple 
to and relax qubits. They considered three models: (I) 
a single electron trap in tunnel contact with a metallic 
gate, (II) a single electron occupying a double trap, and 
(III) a double trap that can absorb/emit a Cooper pair 
from the qubit or a superconducting gate. In all models 
a uniform distribution of the trap energy levels was as- 
sumed. One, then, can show that the distribution for the 
two-level systems corresponding to the models II and III 
are linear in the energy level splitting, P(e) oc e. Since 
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the switching in these models is tunneling dominated, we 
find that -P(A) oc 1/A. Therefore, both models II and 
III are characterized by distribution P(e, A) = Ac/ A, de- 
scribed above, and hence can naturally account for the 
experimentally observed low- and high-frequency noises. 
In contrast, in the model of uniformly distributed single- 
electron traps (model I), we find that for small tunnel 
rates, the high-frequency noise is inversely proportional 
to frequency rather than Ohmic Q. 

In this letter we did not address the question of the 
statistics of the low-frequency noise, nor the associated 
problem of a particular decay law of the dephasing pro- 
cess. These statistics will depend on the distribution of 
the coupling strengths Vj. For certain distributions the 
individual strongly coupled fluctuators may be impor- 
tant [2^, [2^, and the statistics is non-Gaussian. For 
ensembles of many weakly coupled fluctuators Gaussian 
statistics emerges 

To conclude, we have shown that an ensemble of co- 
herent two-level systems with the distribution function, 
P(e, A) oc e/A, produces Ohmic high-frequency noise 
and, at the same time, 1/f low-frequency noise with 
strength which grows with temperature as T 2 . The two 
branches of the noise power cross at uj ~ T in accor- 
dance with the experimental observation A relation 
between low- and high-frequency parts of the spectrum is 
a general property of ensembles of coherent systems with 
discrete energy levels. 
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